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INTRODUCTION 

Local, or short-range, order in disordered alloys is an important and exciting phe- 
nomenon which is quantified in electron, X-ray and neutron scattering experiments. 
It is discussed in many excellent reviews and books, 1 " 6 as well as in the multitude of 
original research papers. 

This relatively short review of the subject does not attempt to discuss all aspects of 
the problem of local correlations in alloys. In particular, we will not touch such issues 
as multiatom (cluster) interactions, static displacements and vibrations of alloy atoms, 
partially ordered, multicomponent or amorphous alloys. As a result, we will concentrate 
on the Hamiltonian traditional for the considered problem, that of the Ising model on 
a rigid ideal lattice with pair, but otherwise arbitrary (i.e., of any range) interatomic 
interactions. 

The central object of the paper is the pair correlation function of the corresponding 
dynamical variables of the model, the occupation numbers or spin variables, the Fourier 
transform of which is proportional to the intensity of diffuse scattering caused by atomic 
short-range order. The main aim is to show that the expression for this quantity has 
certain internal structure analogous, e.g., to that of the averaged Green's function 
used in the electronic theory of disordered alloys. This structure is independent of the 
approximation used for the quantitative description of correlations. As will be seen, 
this structure alone, without further specification of a particular theory of short-range 
order, allows us to see new possibilities in diffuse scattering, some of which have recently 
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been observed experimentally. 

The present contribution is organized as follows. First two sections are auxiliary 
and serve mainly for reference purposes; the former introduces relation between the 
phenomenological Hamiltonian of binary solid solutions and the Ising model, as well 
as necessary definitions and formulae, while the latter describes briefly standard ap- 
proaches to the theory of electronic structure of disordered alloys. Readers familiar 
with the material contained in these sections can skip it and proceed directly to the 
third section where the key expression ( |75|) for the diffuse intensity is introduced. In 
the fourth section its derivation and relation with alloy thermodynamics are discussed. 
Without any doubt, many readers would quickly realize that Eq. ( |ToD is simply one 
of the possible forms of the famous Dyson equation. Such readers can then focus on 
the two last sections. The fifth section reviews, from the point of view adopted in this 
paper, existing theoretical approaches - both traditional and relatively new - to the 
problem of calculation of the pair correlation function and the diffuse-scattering inten- 
sity. Finally, in the sixth section based on the recent author's work it is shown how 
the internal structure of Eq. ( j75|) leads to understanding existing and predicting new 
diffuse-scattering anomalies. 



LIST OF ABBREVIATIONS 



AE 


alpha-expansion 


ATA 


average t-matrix approximation 


CPA 


coherent potential approximation 


GEM 


gamma-expansion method 


HTE 


high-temperature expansion 


KCM 


Krivoglaz-Clapp-Moss (approximation) 


MC 


Monte Carlo (simulation) 


MFA 


mean-field approximation 


PCF 


pair correlation function 


RPA 


random-phase approximation 


SM 


spherical model 


SRO 


short-range order 


SSA 


single-site approximation 


VCA 


virtual crystal approximation 



DESCRIPTION OF A BINARY ALLOY 

We will consider the standard model of a binary alloy used in the statistical theory 
of ordering. 1-6 In this model two sorts of atoms (A and B) are distributed over N sites 
of a rigid lattice; there are no vacancies or other lattice defects. For simplicity reasons 
the consideration is confined to the case of the lattice with one site per unit cell. All 
sites of the lattice are equivalent, i.e., there are no distinct sublattices; this situation 
corresponds to disorder or ferromagnetic ordering. A particular alloy configuration is 
fully described by the set of occupation numbers p", 

a _ J 1, atom of type a at the lattice site i, , > 

^ l [0, otherwise. 
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The occupation numbers satisfy obvious relations 



= 1, (2) 

a 

E^r = at-, (3) 

i 

Y^N a = N , (4) 

where iV a is the total number of atoms of type a, and the last equation holds because 
the total number of alloy atoms is equal to the total number of lattice sites. Indices a 
and (3 take only two values, A or B. It is assumed that the interatomic interactions are 
pairwise; an atom of type a at site i interacts with an atom of type f3 at site j with a 
potential V°f . Then the configurational part of the alloy Hamiltonian is 

ff^EEO^' ( 5 ) 

ij a/3 

The main aim of statistical mechanics is the calculation of the partition function 7 

Z = Trexp(-[3H) , (6) 

where (3 = l/k^T, ks is the Boltzmann constant, T the absolute temperature, and 
Tr denotes the trace of a matrix. It is convenient to avoid restriction (0) on the total 
number of atoms of each type and work in the grand-canonical ensemble, calculating 
the grand partition function Z, 

where fi a is the chemical potential of atoms a. In this approach the number of atoms is a 
function of the corresponding chemical potential, and after having done the calculation 
the chemical potential is adjusted to get the required number of atoms. To unify 
canonical and grand-canonical treatments, i.e., to get the same formula (||) for Z, an 
effective Hamiltonian H is often introduced: 

H = H-J2^ a N a . (8) 

a 

The expression for the effective Hamiltonian is then 

# = ^EE vfptv^ -EE « • (9) 

ij a/3 i ol 

As is well-known, the statistical-mechanical problem defined by the Hamiltonian 
([D can be transformed into the equivalent problem of an Ising magnet in a magnetic 
field. The occupation numbers are not yet all independent, due to Eq. (Q). Dependent 
variables can be eliminated by introducing a spin variable Sj, 

J +1, spin up (atom A) at the site i, , , 

\ — Ij spin down (atom B) at the site i. 

The relations between the occupation numbers and the spin variables are 

Pt = \(l + Si), (11) 
Pf = (12) 



Inserting Eqs. (|TTJ) and ( |T2| ) into the expression ([5]) for the effective Hamiltonian H, we 
find that, apart from the configuration-independent term, the Hamiltonian is that of 
the Ising model, 

H = --J^JijSiSj - hJ^Si , (13) 

ij i 

where the effective exchange integral and the magnetic field h are given by 

Jij = ~\v is , (14) 
h = . (15) 

^ = l(V4 A + V* B )-V4 B , (16) 

and is the pair ordering potential. The second term in Eq. (|15| ) does not, in fact, 
depend on the index i because all interatomic potentials are functions of only the 
distance between interacting atoms. 

For the following we define concentrations c a and magnetization m which are statis- 
tical averages of the occupation numbers ([]]) and the spin variables (0), respectively, 

c a = <p?>, (17) 
m = ( Si ) , (18) 

getting from Eqs. ( |TT] ) and (|T 



c=i(l+m), (19) 

where c = c A = 1 — c B . Here the symbol (...) denotes the statistical averaging with 
the effective Hamiltonian H, Eq. (0) or ([IB]), 



Tr 



exp 



(-PH) 



Trexp(-/3 J ff) ' ^ 
We also introduce irreducible PCFs for the occupation numbers, 

G% P =(P>I)-(p-)($) , (21) 

and for the spin variables, 

Gij = (siSj) - (si) (sj) . (22) 

The notion of irreducibility comes from the diagram technique for the Ising model, in 
the framework of which it can be shown that the irreducible PCF does not contain 
so-called split diagrams which are present in the expansion for the correlation function 
(siSj). 8 Using relations ( |TT| ) and ([12]) between the occupation numbers and the spin 
variables, it is easy to obtain the correspondence between the PCFs (pl|) and (p^), 

riAA _ r~iBB _ r^AB _ n /oq \ 

and, therefore, there exists only one independent PCF for the occupation numbers. In 
the theory of alloys the Warren-Cowley parameters a^, instead of the PCFs (pT]), are 
traditionally used for the description of SRO. They are defined as 



a ., = 1 _ (vtvf) m 

a t3 c (i- c )' U1) 
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and, according to Eqs. (0), (|IJ), (|2T|) and (|23|), 



Gij = 4c(l — c)ctij = (1 — m )atij . (25) 



The diagonal matrix elements of the PCFs (^1|) and ( [22D can be expressed in terms 
of the averages (|17D and flT8|). Taking into account identities (p?) 2 = and sj = 1 



which follow from the definitions (|l|) and (fLOf ), we have 

= GS fl = -G^ = c(l - c) , (26) 
= 1 - m 2 , (27) 
a« = 1 . (28) 

The last equation leads to the important sum rule in the reciprocal space. Let us 
introduce the lattice Fourier transformation, 

/(k) = ^£^exp(*kr 4 ,) . (29) 

ij 

fij = ^J dk /( k ) eX P (~ ikr ij) > ( 30 ) 

where is an arbitrary lattice function, k the wavevector, = r,, — r^, r, the radius- 
vector of the site i, and the integration in Eq. ( p0|) is carried out over the Brillouin zone 
of volume Q. Eq. (|28|) can then be rewritten as 



1 



dk a(k) = 1 . (31) 



Recalling that the Fourier transform a(k) of the SRO parameters a^- is the SRO part 
of the diffuse-scattering intensity measured in Laue units, 5 we conclude that the sum 
rule fl3T|) expresses the property of conservation of the integrated SRO intensity. 



ELECTRONIC THEORY OF DISORDERED ALLOYS 

In what follows we will occasionally turn to ideas which were developed in the 
theory of electronic structure of disordered alloys. 2 ' 6,9-11 With this in mind, we deviate 
now from the main theme of this paper and consider briefly (and rather formally, 
without discussing the physical meaning of derived formulae) methods of calculating 
the Green's function of electrons averaged over possible alloy configurations. There is 
some overlapping of notations used in this section with those in the rest of the paper; 
however, the present section is quite isolated from the statistical-mechanical part of the 
discussion, and hence this overlapping should not lead to any confusion. 

The central role in the one-electron theory of disordered alloys is played by the elec- 
tronic Green's function G which may be defined as a resolvent of the alloy Hamiltonian 
H, 

G=(E-H)- 1 , (32) 

where E is the energy of an electron. The Hamiltonian of a disordered binary alloy is 
usually assumed to be a sum of two operators, 

H = H + V . (33) 

The first of these two terms, Hq, is translationally invariant, while the second, the ran- 
dom one-electron potential V, depends on particular alloy configuration and is usually 



r 



assumed to be diagonal in the site representation. The operator V is a sum of individual 
potentials Vi centred at each site i and acquiring two possible values, V A and V B , in 
accordance with the type of atom occupying the considered site. 

Propagator expansion 

The potential V is often viewed as a perturbation, though not necessarily small, of 
the initial unperturbed Hamiltonian Hq. Defining the unperturbed Green's function, 

G = (E-H )- 1 , (34) 

one can easily construct the Dyson equation, 

G = G + G VG . (35) 

Iterating this equation and averaging over all possible alloy configurations (the averages 
are denoted by brackets), we get the propagator expansion 

(G) = G + G (V)G + G (VG V)G + ... . (36) 

The averaged Green's function satisfies another Dyson equation, 

(G) = G + G J:(G) , (37) 

where the operator £ is called the self-energy. The self-energy is, in terms of the 
diagrammatic expansion for the averaged Green's function (G) generated by Eq. 
the irreducible part of (G), i.e., the sum of all graphs for (G) which cannot be separated 
into two parts by cutting a single bare-propagator line Gq. From Eq. (|37|) it follows 
that 

(G) = (Go" 1 - S)- 1 = (E-H - S)- 1 , (38) 
the second equation comes from Eq. (p 7 !]). 

Locator expansion 

There exists another perturbation series, the locator expansion, which is the expan- 
sion in powers of the unperturbed Hamiltonian H rather than the potential V. In this 
case the unperturbed Green's function is that of the localized atomic states, 

g = (E - V)- 1 , (39) 

and the corresponding Dyson equation has the form 

G = g + gH G. (40) 

As in the case of the propagator expansion, we iterate this equation and average over 
the ensemble of configurations term by term, obtaining the following expansion: 

(G) = (g) + (gH g) + (gH gH g) + ... . (41) 

By introducing the so-called fully renormalized interactor U, 

U = H + H GH , (42) 
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full formal analogy between the propagator (Eq. fl3"6|)) and locator expansions is achieved: 

(U) = H + H (g)H + H (gH g)H + ... . (43) 

The corresponding sum of all irreducible graphs a in the case of the locator expansion 
is called the locator; similarly to Eq. fl3"T|), we have 

(U) =H + H a(U) . (44) 

From Eqs. (|42|) and ([44]) it immediately follows that in terms of the locator the Dyson 
equation and the expression for (G) are 

(G) = a + aH (G) , (45) 
(G) = (a-'-Ho)- 1 , (46) 

respectively. The relationship between the locator a and the self-energy S, according 
to Eqs. © and fl4§), is 

a = {E-Y,y l . (47) 

Renormalization of one-electron potential 

We now renormalize the potential V subtracting a configuration- independent, site- 
diagonal operator S and adding it to the unperturbed Hamiltonian H : 

H = H + V = (H + S) + (V-S) = H + V . (48) 

If one defines the unperturbed Green's function with respect to the renormalized Hamil- 
tonian H, 

G=(E- H)- 1 , (49) 

then the Green's functions G and G are related by the Dyson equation analogous to 
Eq. ©, 

G = G + GVG . (50) 



As follows from Eqs. (g8|) and (|48|) 



(G) = (E-H-t)- 1 , (51) 

where E relates to H and V in the same way as S to H Q and V, and 

£ = S + £. (52) 

The operator S may be regarded as an initial approximation for the exact self-energy E. 

The next step is to introduce the total scattering operator T and express the exact 
Green's function G and the self-energy S in terms of T. The total scattering operator 
is defined by the relation 

T = V + VGT , (53) 

which gives 

T = (V- 1 - G)- 1 . (54) 

Excluding V from Eqs. ([5(]) and ([54]) and averaging the result over realizations of the 
random potential, we get the relation between (G) and (T): 

(G) = G + G(T)G . (55) 
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Taking into account the relation (^) between the averaged Green's function and the 
self-energy and using Eq. (|55|), we finally obtain 

Z = S + ((T)- 1 + G)- 1 . (56) 

Thus, having calculated the averaged total scattering operator (T) one can determine 
the averaged Green's function (G) and the self-energy £ according to Eqs. (|55| ) and 
(|56|). The problem of the description of the disordered alloy is therefore reduced to the 
problem of finding reasonable approximate expression for the operator (T) . The opera- 
tor T can be expressed in terms of operators describing scattering on individual atomic 
potentials V{. To do this, it is convenient to decompose the Green's function G into 
two parts which are diagonal and off-diagonal in the site representation, respectively: 

G = Gd + G d ■ (57) 



Inserting this equation into Eq. (|54|), we have after some straightforward algebra, 

T = (r 1 - Go*)- 1 , (58) 

where the operator 

t = (V- l -G d y 1 (59) 

is, similarly to the potential V, site-diagonal and represents a sum of individual scat- 
tering operators t% corresponding to atomic potentials. 



Single-site approximations 



At any level of approximation there exist two main approaches to the problem of 
calculation of the averaged total scattering operator (T) and, consequently, the averaged 
Green's function (G) and the self-energy S. First, it is possible to choose the operator S 
from the very beginning, and then calculate these three quantities using Eqs. (p^)-(p5|); 
this is the non-self-consistent scheme. In most cases the VCA choice S = (V) is used. 
Second, one can consider S as an operator variable and, noticing that the scattering 
operator T is a function of S, determine the latter as a solution of the equation 

(T[S}) = . (60) 

This equation gives, according to Eqs. Q55|) and (|56|), 

(G) = G, (61) 
£ = S , (62) 

and such an approach is called self-consistent. To make the calculation of the av- 
eraged total scattering operator practically possible, the following decoupling of the 
configurationally-averaged Eq. (pB[), called the SSA, is usually adopted: 

(T) = ((r 1 - G od y l ) S -H ((t)- 1 - G^y 1 . (63) 

In the framework of the SSA Eq. (|56D for the self-energy takes the form 

S^ = 5 + ((t)" 1 + G (i )- 1 , (64) 

which means that in this approximation the self-energy is diagonal in the site repre- 
sentation. Consider now the two approaches in combination with the SSA. The non- 
self-consistent one, as is seen from Eqs. ([59]) and (|64"D , gives the self-energy as a sum of 



G 



the operator S and the effective scattering potential V e ff corresponding to the average 
scattering operator (t), 

Veff = ((tr 1 + G d y l ; (65) 

accordingly, this method of calculation is known as the ATA. The application of the SSA 
to the self-consistent scheme (Eq. (|60|)) gives the following equation for the evaluation 
of the self-energy: 

(t[S]) = 0. (66) 

This equation determines the CPA. 



EXPRESSION FOR SHORT-RANGE ORDER DIFFUSE INTENSITY 

We now derive a formally exact expression for the SRO part ct(k) of the diffuse- 
scattering intensity. A key observation here is that the spin PCF ( ^2|) satisfies the 
Dyson equation, 8 

G = a + a AG , (67) 

where A = (3 J and, like in the previous section, matrix notations are used. Here a 
is again the sum of all irreducible graphs in the diagrammatic expansion for G, but 
irreducibility is now defined with respect to the interaction-to-temperature ratio A. 
Eq. flBTD has the same form as the Dyson equation fl4"5] ) for the Green's function of 
electrons (|32|) averaged over alloy configurations. From Eq. (^) it follows that 



G = [a' 1 - Aj . (68) 
This expression for the spin PCF is analogous to Eq. ([46]). 



The irreducible part a of the PCF G is sometimes called the self-energy. 8 However, 
to maintain the analogy with the electronic theory of alloys (i.e., with the terminology of 
the previous section) we will refer to this quantity as the locator, and reserve this term 
for another object, defining the PCF self-energy £ by the relation similar to Eq. fl4"T|): 

a = -5T 1 . (69) 

In terms of S Eq. ( |68"D becomes 

G = (— £ - A)" 1 . (70) 

We also introduce the locator a and the self-energy £ for the occupation-number PCFs 
G«*(Eq. (0)): 

* = (71) 
£ = 4S . (72) 

Then matrices a and S are related by the same Eq. 

a = -XT 1 , (73) 

and we have 

Q AA = qBB = _ Q AB = c(1 _ c)a = + 2/3V)- 1 . (74) 



n 



Written in k-representation, this equation leads to the following expression for the 
diffuse-scattering intensity: 12,13 

a(k) = T . (75) 



c(l-c) 



-E(k) + 2pV(k) 



The central quantity of interest for us here is the PCF self-energy E(k); apart from 
the wavevector, it depends also on two other variables, temperature and concentration. 
Later on we will be focusing on existing approximations for the PCF self-energy. 



SELF-ENERGY AND THERMODYNAMICS 

In the last section the Dyson equation ( p7|) and the related expression ( [70]) for the 
spin PCF were simply postulated. However, it would be useful to know how Eq. ( |70| ) 
could be derived and how the self-energy is related to the thermodynamics of the sys- 
tem. 

Variational formulation of statistical mechanics 

As was mentioned before, the main task of the statistical-mechanical treatment is to 
calculate the partition function (|6]) of a system described by the Hamiltonian H. The 
Hamiltonian usually is a linear combination, 

H = J2x n a n , (76) 

n 

of some operators a n with coefficients x n . Variables a n conjugated to the parameters 
Xji are defined as averages (^) of the operators a n , 

dF 

o: n = (a n ) = t; — , (77) 

where 

F = -k B T\nZ (78) 

is the free energy of the system, and the second equation in ( |77D comes from the 
definition (|) and Eq. (|76|). Our aim now is the calculation of the averages a n and the 
free energy F as functions of the parameters x n . This problem can be formulated as a 
variational one, if the Legendre transform T(a) of the free energy F(x) is introduced: 14,15 

T(a) = F(x(a)) - ^a n x n (a) . (79) 

n 

Here the averages a n are the independent variables, and x n (a) are solutions of Eqs. ([FT]). 
Differentiating T(a) with respect to a n and using Eqs. ([H]), we get 

^ = Y, — — -Y,(— X rn + a m ^) =-X n . (80) 

dotn m 0x m dot n m y dot n d(y n J 

Finally, introducing another function $(x,a), 

<f>(x,a) = r(a) + ^2<y n x n , (81) 



i n 



which depends on both a n and x n , we find from Eq. fl8"0|) that it is stationary with 
respect to variations of a n at fixed x n : 

. 82 



At the stationary point a n = a n (x), where a n (x) are solutions of Eqs. (|S0|), a) 
coincides with the free energy F(x): 

a(x)) = F(x) . (83) 

Function $(x, a) is usually called the variational free energy. Noting that the internal 
energy is the statistical average fl20|) of the Hamiltonian, E = (H), &(x, a) can be 
written in the standard thermodynamic form: 

$(x,a) = E(x, a) - TS(a) , (84) 

E(x,a) = ^x n a n , (85) 

n 

5(a) = -(3T(a) , (86) 

where E(x, a) and 5(a) are the variational internal energy and configurational entropy, 
respectively. Eq. (|82|) now becomes 

T^ = x n . (87) 



da 



n 



From Eqs. (tfTf) and (g) it follows that 

d 2 T d 2 F _ 3x n dai _ dx n _ ^ 
i da n dai dxidx m ^ da t dx m dx m 

Eq. ( ^8|) shows that matrices of second derivatives of the free energy F(x) and its Leg- 
endre transform T(a) (or the variational free energy $(x, a) which differs from T(a) 
only by the bilinear term Y^n x nOt n ) are mutually inverse up to a sign. 

First Legendre transformation for the Ising model 

We will now apply the general technique of the Legendre transformations outlined 
above to the particular case of the Ising model in the inhomogeneous magnetic field. 
The Hamiltonian of the model is a straightforward generalization of Eq. (|T3|): 



H = -- J2 JijSiSj - hiSi . (89) 

ij i 

From comparison of Eqs. (|76|) and ( |8~9~D it follows that the latter contains two kinds of 
operators a n - single spin variables Sj and products SiSj of two spin variables. Corre- 
sponding parameters x n , except for sign, are hi and Jy, and variables a n conjugated 
to these parameters are (sj) = m; (Eq. (|T8|) ) and (sjSj) = m^mj + (Eq. (p^)). In 
the general case considered before the Legendre transformation was carried out with 
respect to all parameters result, all conjugated variables a n were calculated 

using the variational principle. Here, however, this approach will be applied only to 
the magnetic field hi, and only the magnetization m; will be calculated by means of the 
variational procedure. The resulting partial transformation is called the first Legendre 
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transformation. 15 17 In this case hi plays the role of x n , and comparison with Eq. ( |76| ) 
shows that a n corresponds to — Sj. Eqs. (|77|), (|79"1)-(|83"D and (|8|) now become 

r(m,J) = F{h(m),J) + ^mihi{m) , (90) 

i 

$(/i,m;J) = r(m, J) -^mA , (91) 

i 

-m, , (92) 

hi, (93) 

= °- < 94 » 

<$>(h,m(h);J) = F(h,J) , (95) 
y d 2 T d 2 F = <9 2 $ d 2 F _ 5 



dF 

dhi 
dY 



where, in the same manner as earlier, hi{m) and rrii(h) are solutions of Eqs. (|9~2"D and 
(j93[), respectively. 

Derivation of the Dyson equation and meaning of the self-energy 

First of all, we note that from the definitions of the partition function, the statistical 
average, the PCF and the free energy (Eqs. (^), (|2"0"D, ( F2"2"D and (|T8D) it follows that for 
the Hamiltonian (|89"D 

d 2 F 

aim- = ' < 97) 

Then, combining this equation with Eq. fl9~6|), we obtain 



G « = fcsT tt-jt- • ( 98 ) 

\omom J „ 

This notation means that the real-space matrix element of the PCF is proportional to 
the corresponding matrix element of the inverse of the matrix whose matrix elements 
are second derivatives <9 2 $ / 'dmidrrij of the variational free energy with respect to the 
corresponding magnetizations. 

The variational free energy can always be written as a sum of its mean-field (Bragg- 
Williams) part and the non-mean-field correction: 

$ = <5> MFA + 5$ , (99) 

qMFA = E MFA _ TS MFA _ ^qq) 

The expressions for the mean-field internal energy E MFA and the configurational en- 
tropy S MFA are well-known: 1-6 

e mfa = ~J2 J 'i m " n i-T, h i m < . < 101 ) 

= -feWi^lnl^ + l^ilni^i) . (102) 



Substituting Eqs. (ggJ-QIDg) into Eq. @ and noticing that 



d 2 E MFA j 
dmidrrij 13 

we recover the result ([70]) for the PCF (and, therefore, the Dyson equation ([H])) with 
the following expression for the self-energy: 



gigMFA Q2 (gfy 

dmidrrij dmidrrij 



>:, ; = if Q „_ - /? L 1„ • (ioi) 



It is seen that the self-energy is the sum of the second derivatives, with respect to the 
magnetizations, of the two terms contributing to the expression for the variational free 
energy $: the mean-field configurational entropy and the non-mean-field part of $. 
Noting further that 



d 2 S MFA g 



dmidmj 1 — mf ' 



(105) 



we finally obtain 



Ei . = __^ p 92 W . (106) 

13 1 — mf dmidmj 

The first term in this expression is diagonal in the direct space. This means that, 
back to the homogeneous case in which all the lattice sites are equivalent, this term is 
k-independent in the reciprocal space, and all the wavevector dependence of the self- 
energy comes from the second term, i.e., is the result of the non-mean-field corrections 
to the MFA variational free energy. 



APPROXIMATIONS FOR THE SELF-ENERGY 

In this section our attention will be focused on the formally exact result ( [75] ) for the 
SRO diffuse intensity a(k), which is the Fourier transform of the Warren-Cowley SRO 
parameters ajjj (24). Available theories of SRO will now be considered, in the light of 



the structure of Eq. ( |7q) , as different approximations for the self-energy S. 

Random-phase (Krivoglaz-Clapp-Moss) approximation 

The simplest and by far the most popular theory of SRO is the RPA (or, in the 
alloy language, the KCM approximation) : 5,18,19 

aRPA ^ - 1 + 2c(1 ! c)mk) • m 

Eq. (p.07|) is usually derived using mean-field-like arguments. Comparing Eqs. fl75D and 
(|1 07]) , we see that 

£ "< k » = -^7) ■ < 108 > 

The RPA self-energy ( |108|) is thus wavevector- and temperature-independent; it de- 
pends only on alloy composition. Returning via Eqs. (19) and (|T2|) to the magnetic 



language used in the previous section, we conclude that the RPA for the self-energy 
corresponds precisely to neglecting the second term in the right-hand side of Eq. ( |106| ). 
The non-mean-field contribution to the self-energy is therefore ignored in the RPA, 
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and the self-energy in the direct space is simply the second derivative of the mean-field 
configurational entropy with respect to the magnetization: 



^ Q2 gMFA 

^ drriidrnj 1 — mf ^ 

From the point of view of the terminology used in the electronic theory of disordered 
alloys, the RPA for the self-energy could be referred to as the SSA; indeed, both the RPA 
(Eq. (|104|) ) and the SSA (Eq. QSij) ) self-energies are diagonal in the site representation. 
The RPA resembles the non-self-consistent SSA, i.e., the ATA, in the sense that both 
approximations define the corresponding self-energies explicitly. 

The RPA reduces to the well-known Ornstein-Zernike description of correlations, 7 
when only those wavevectors which are close to the position ko of the absolute minimum 
of the interaction V(k) are considered (this approximation corresponds to the long-wave 
limit in the case of ferromagnetic ordering). Let us expand V(k) in powers of q = k — ko 
and retain only the lowest-order (quadratic) term, 

V(k)^V(k ) + W9ijQiQj , (HO) 

where g is the 3 x 3-matrix of second derivatives of V(k) at k . We will take only the 
simplest example of cubic symmetry, when = g5ij] in this case J2ij dijQiQj — 9Q 2 -, 
where q = |q|. Substituting the result into Eq. ( |107| ), we obtain 

aRPA M = z (T T ff n x i , (in) 

k B {T - T c ) + c(l - c)gq 2 

where it is taken into account that the mean-field result for the instability temperature 
T c is 

T c = 2c(l-c)\V(k )\/k B ; (112) 

at the position of the absolute minimum the interaction value is negative, since the 
average of the interaction over the Brillouin zone is zero (see Eq. ( |116|) below). In real 
space we get asymptotically (i.e., at large distances) 

I 

aff A oc — exp(-ik r ij - r^/^) , (113) 
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'j 



where = |i*y|, and the correlation length ^ is 



C(1 - 0)9 (114) 



\ k B (T-T c 



Eqs. ( |111| )- (|114[) represent the Ornstein-Zernike result for the PCF. 



However, the RPA expression ( |107| ) has a serious disadvantage: it is unable to satisfy 
the sum rule fl3"ID . Using the identity 

1 4c 2 fl - c) 2 B 2 V 2 (k) 



1 + 2c(l - c)/3V(k) K ,y K ' 1 + 2c(l - c)(3V(k) ' 



it can be shown that the RPA formula ( 1Q7|) always leads to the overestimation of the 



value of the integral over the Brillouin zone (|3"T|). Since 

^ J dk V(k) =Va = 0, (116) 



from Eq. ( |115| ) it follows that 



q J [k) ~ 1 + n 7 dk i + 2 C (i - c)^y(k) ' lii7J 

and, therefore, 

i 1 dk a flPA (k) > 1 . (118) 

This integral is close to unity only at sufficiently high temperatures. As temperature 
decreases, the deviation from the value prescribed by the sum rule becomes more and 
more significant, and the integral finally diverges at the instability point. 20,21 



Spherical model 

Another analytical expression for a(k) is given by the SM, 6,19,22 ~ 24 also known as 
the Onsager cavity field theory, 25,26 

, SM ' 



a bM (k) = r , (119) 



c 1 



•E SM + 2/3V(k) 



where E is, at fixed temperature and concentration, a number determined from the 



sum rule (|31|). Therefore, the sum rule is satisfied in the SM by definition, contrary 
to the case of the RPA. From the definition of the SM it also follows that the self- 
energy depends not only on concentration, like its RPA counterpart ( |108| ), but also on 
temperature. Nevertheless, the SM self-energy is still wavevector-independent. The 
explicit expression for S 5M can be derived from the sum rule (fBlf): 



t SM = t RPA + 5t , (120) 
5t = 2l3'£a ij V ij = ^Jdka(k)V(k). (121) 

3 

Similarly to the RPA, the SM is analogous to the SSA, since the SM self-energy is 
diagonal in the direct space. However, the SM is rather the self-consistent SSA, like the 
CPA, because the sum rule here plays the role of the CPA self- consistency condition. In 
fact, the sum rule is the self-consistency condition: it simply means that the diagonal 
matrix elements of the approximate and exact PCFs are the same. More surprisingly, it 
was shown 27 ' 28 that both the CPA and the SM could be obtained by summation of the 
same sets of diagrams in the corresponding perturbation expansions. We can conclude, 
therefore, that the SM is the CPA for the Ising model. 



Cluster variation method 



The CVM is at present the standard technique for quantitative calculation of 
thermodynamic properties of alloys. It is discussed in great detail in almost every book 
or review on the subject. 1-3 ' 6 We do not attempt to do this here; instead, we will 
consider only those features of the CVM which are relevant to our discussion, without 
going into technical aspects of the method. 

The CVM is essentially a procedure which allows us to derive an approximate ex- 
pression for the variational configurational entropy S(a) of the system. The CVM 
entropy is a function of probabilities of various atomic configurations on lattice clusters 
which belong to the so-called basic, or maximal, cluster. A particular CVM approxi- 
mation for S(a) is therefore defined by the choice of the basic cluster. Combined with 



i r-. 



the variational internal energy E(x, a), the CVM configurational entropy gives the ex- 



pression (84) for the variational free energy <S>(x,a). The operators a n in the alloy 
Hamiltonian are products of the spin variables or of the occupation numbers, and the 
averages a n (Eq. ( [77|) ) are thus related to, or coincide with, cluster probabilities enter- 
ing the expression for the CVM configurational entropy. The variational free energy $ 
is then minimized with respect to all cluster probabilities, taking into account various 
self-consistency constraints. 

The self-energy obtained in the framework of the CVM depends, in general, on all 
three parameters - temperature, concentration and wavevector. The problem with the 
CVM, however, is that this method is intrinsically numerical and does not lead to an- 
alytical approximations for pair correlations. The reason for this is that in most cases 
the number of cluster variables used to get a reasonably accurate formula for the con- 
figurational entropy is far too large. From the point of view of the general technique of 
the Legendre transformations, the CVM corresponds to the high-order transformation 
with respect to all coefficients x n conjugated to cluster variables which are involved in 
the expression for the variational configurational entropy S(a). Therefore, as far as 
SRO is concerned, general Eq. fl5S|) is still valid, as is the Ising model-specific Eq. (|9^|). 
However, in practical sense this case is very different from that of the first Legendre 
transformation. In the latter, the inversion of the matrix <9 2( 3> / 'dada is carried out triv- 
ially using the Fourier transformation. In the CVM this object in the reciprocal space 
is still a sufficiently large matrix, and in all realistic situations it is necessary to resort 
to numerics. Analytical formulae were obtained only for such simple model cases as the 
pair (also known as quasichemical, or Fowler-Guggenheim, or Bethe) approximation, or 
the square (Kramers- Wannier) approximation for the nearest-neighbour Ising model on 
the square lattice. 2 ' 20 Besides, the CVM suffers from the same drawback as the RPA, 
though to a lesser extent: the integrated intensity is not conserved, and its behavior 
with temperature is similar to the case of the RPA. 2 ' 20 

High-temperature expansion 

As we have seen, the RPA and the SM are the analogs of the SSAs in the electronic 
theory of disordered alloys and thus fail to take account of the wavevector dependence 
of the self-energy. On the other hand, the CVM leads to the k-dependent self-energy, 
but loses the simplicity of the former two approximations by not providing analytical 
expressions for correlations. At the semiquantitative level of the RPA this problem can 
be cured by using the HTE for the self-energy. 30 

To do this, we return to Eq. ( |104j ) which, in combination with Eq. ( |109| ), can be 
written as 

£ = Z RPA + <5£ , (122) 
d 2 (5$) 

*E« = -P^jr 1 - (123) 
avrLiOmj 

It is known that the RPA is exact to first order in 1/T (i.e., in A); 2 it means that 
5$ and <5£ are of order A 2 . It is also known how to construct the first Legendre 
transformation T (and, therefore, 5$) for the Ising model. 15-17 Sorting all contributions 
to 5$ (diagrams) according to the number of lines A, we obtain the HTE; first eight 
orders in A are available in the literature. 17 For simplicity reasons, the discussion here 
is confined to two first orders, 

-^4° + T2^ + l A +0(A 4 ), (124) 
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where a line corresponds to A, and a vertex with n attached lines represents the function 
u n (rrii). For the diagrams in Eq. ( |124|) n is equal to either 2 or 3: 

u 2 {mi) = 1 - m\ , (125) 
u 3 K) = -2m*(l-m?) . (126) 

The expressions for the diagonal and off-diagonal parts of ST thus are, according to 
Eq. O, 



<5£ 



104l^v +oi4,) - < i27) 

= \ Oi + g t©t + /S, + °< a4 ) • ^j- ( 128 > 

Here a vertex with n internal lines and k external legs corresponds to the kth derivative 
of the function u n (m). The corresponding analytical expressions have the form 

SEa = -£(l-m?)A£- 477N£m,(l-m?)A3i 
i i 
-£(l-m 2 )(l-m 2 )A jfe A^ + 0(A 4 ) , (129) 

kl 

ST^ = 2m,m,A^. + ^(l-3m J 2 )(l-3m J 2 )A| 

+4771*77^ ^(1 - mf)AuAji + 0(A 4 ) , z ^ j . (130) 
i 

In the homogeneous case m; = m for all sites z. Defining constants 

Qi = ^ J a , 0-2 = X! J a > °3 — ^JikJuJki (131) 

and lattice functions 

(A)ij = ^ij , C/^ij = ^ , (/3)y = Jij^JilJjl , (132) 

which depend only on the interaction J, one can finally write the expression for ST in 
the k-space: 

<5£(k) = -(1 - m 2 )f3 2 ai + 2m 2 /3 2 /i(k) -4m 2 (l - m 2 )/3 3 a 2 
-(l-m 2 ) 2 /5 3 a 3 + ^(l-3m 2 ) 2 /3 3 / 2 (k) 

+4m 2 (l - m 2 )/3 3 / 3 (k) + 0(/3 4 ) . (133) 
In alloy notations this result reads 



ST{k) = -4c(l-c)/? 2 ai + 2(l-2c) 2 /? 2 /i(k) + 8c(l-c)(l-2c) 2 /3 3 a 2 
+8c 2 (l - c) 2 /3 3 a 3 - |[1 - 6c(l - c)] 2 /3 3 / 2 (k) 

_ 8c (l - c )(l - 2c) 2 /3 3 / 3 (k) + 0(/3 4 ) , (134) 



where 



a! = > «2 = E^?' a3 = Y, v * v iMa> ( 135 ) 

(fih = \'0- (hh = v 5> (fa)v = V ij y EVuV jl . (136) 



Alpha- and gamma-expansions 



The approximate expressions (|133|) , (|134j) for the self-energy given by the first orders 



of the HTE satisfy almost all requirements: they are analytical, relatively simple and 
take into account dependence on all relevant parameters, including the wavevector. 
However, the limits of applicability of the HTE are essentially the same as those of the 
RPA; the HTE is quantitatively correct only at reasonably high temperatures, as can 
be concluded already from the name of the expansion. What is needed, therefore, is 
an approximation which would combine all the advantages of the HTE self-energy with 
the applicability at moderate temperatures, including the range not far away from the 
transition or instability points. 

The theory of SRO which will now be discussed 12,13 is based on the fairly gen- 
eral procedure 27 of self-consistent renormalization of the bare propagator A -1 in the 
functional-integral representation of the generating functional for correlation functions. 
The resulting expansion for the matrix elements of the self-energy is in powers of the 
matrix elements of the fully dressed propagator, which in the case of the Ising model 
coincides with the PCF Q70|). Two first non-zero orders of this expansion for the off- 
diagonal part of the self-energy were calculated, 12 ' 13 



E y = aG% + bG% + 0{G A ) , i^j, (137) 

(138) 



2m 2 



(1 -m 2 ) 4 ' 
2[(1 -3m 2 ) 2 - 12m 4 ] 



3(1 — m 



2\e 



(139) 



in terms of the alloy variables 



Ey = ~aa% + ba% + 0(a 4 ) , i ^ j , (140) 

\2 



<' = .v ' , 2 t 2 . (141) 



2[c(l-c)P 
[l-6c(l-c)] 2 -3(l-2c) 4 
6[c(l-c)] 3 



(142) 



Eq. ( |140[ ) is the expansion in powers of the SRO parameters, and is therefore referred 
to as the AE, 31 though initially it was obtained in the framework of the GEM. 12,13 
The difference between the AE and GEM is discussed below. In the two calculated 
orders the AE preserves the sum rule ([H|), and the expression for the diagonal part of 
the self-energy then comes from Eqs. (|3l|) and (|140|) , similarly to the case of the SM 
(Eqs. O, O): 

E« = t RPA + 5t u , (143) 
5t tl = 2/?5> ii VS i - £ (aa?. + &4) +0(a 5 ) . (144) 

3 j&i) 



Note that the first term in Eq. (|144|) corresponds to the SM (cf. Eq. (|121 )) which is 



the zero-order approximation for both the AE and the GEM; in the SM the self-energy 
is diagonal (a = b = 0). 

The difference between the GEM and the AE is in the way of selection of leading 
terms in the expansion (|140|) . The GEM, originally proposed by Tokar 27 and further 
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developed by Tokar, Grishchenko and the author, 12 ' 13 ' 32 35 is based on the assumption 
that correlations decrease rapidly with distance, and the GEM expansion parameter is 

7 = exp(-l/0 , (145) 

where £ is the dimensionless correlation length. The leading terms in the diagrammatic 
expansion for the self-energy are selected in the framework of the GEM according to 
the total length of all lines in the diagrams, where the line connecting lattice sites % 
and j represents a^. For example, in the case of three Bravais lattices belonging to the 
cubic system taking into account several first terms of the perturbation theory leads to 
the result 12 ' 13 

t s = aa 2 s + ba 3 s , s = 1 , (146) 
t s = ~ac? s , s = 2,3, (147) 
t s = , s > 3 , (148) 

where subscript s denotes the matrix elements corresponding to the sth coordination 
shell. However, the GEM assumption about the rapid decay of correlations is not always 
valid; e.g., it is incorrect in the cases where distant interactions are essential. The AE 
abandons this assumption and uses instead the SRO parameters themselves as the 
expansion parameters; the leading terms are chosen according to the number of lines 
in the diagrams (i.e., the powers of a^-), since all ay are relatively small. A particular 
AE approximation is defined by neglecting higher-order terms and including only finite 
number of coordination shells in the AE expression ( |140|) for the off-diagonal part of the 



self-energy. The GEM was successfully applied to both the direct and inverse problems 
of alloy diffuse scattering, 12,13,36-38 leading to reliable results everywhere except in the 
vicinity of the instability point, while the AE was used in the analysis of some of the 
diffuse-scattering anomalies discussed in the next section. 31 ' 39 



ANOMALIES IN ALLOY DIFFUSE SCATTERING 

In this last section we will show how expression (|75|) for the intensity ct(k) and, in 
particular, the wavevector dependence of the self-energy S(k) lead to straightforward 
explanation of recently observed unusual features (anomalies) of diffuse scattering from 
disordered alloys and to prediction of some new effects. 

Temperature-dependent Fermi surface-induced peak splitting 

This curious effect (the temperature dependence of the splitting) was discovered in 
1996 by Reichert, Moss and Liang 40 in the first in situ experiment to resolve the fine 
structure of the equilibrium diffuse scattering intensity from the disordered CU3AU alloy. 
The separation of the split maxima changed reversibly, increasing with temperature. 
The same behavior of the splitting was also found 41 by analysing results of the MC 
simulations for the Cu .856Alo.i44 alloy. 42 The fourfold splitting of the intensity peaks 
located at the (110) and equivalent positions (Figure 1) is attributed to the indirect 
interaction of atoms via conduction electrons in an alloy whose Fermi surface has flat 
portions; the effective interatomic pair interaction V(k) itself has split minima in the 
reciprocal space, and their location is determined by the wavevector 2kp spanning these 
flat portions of the Fermi surface. 5 ' 43 It is usually assumed that V(k) is temperature- 
independent. This assumption is justified at least as far as positions of the V(k) minima 
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are concerned, since the 2kp value is unlikely to change over the considered temperature 
range. Besides, the MC calculations 42 in which the increase of the splitting with tem- 
perature was found 41 were carried out for the temperature-independent pair interaction 
parameters. The standard RPA (KCM) treatment (Eq. Ql07| )) predicts that positions 
of the intensity peaks coincide with those of the corresponding minima of V(k), and, 
therefore, the splitting does not depend on temperature. 




Figure 1. 3D reciprocal-space picture of scattering from the FCC alloys discussed in the text. Large 
dots represent the Bragg reflections. Characteristic crosses formed by small dots correspond to the 
split diffuse intensity peaks. 



As will be demonstrated below, this phenomenon can be easily understood by em- 
ploying the notion of the k-dependent self-energy. 31 Let us consider the ct(k) profile 
along one of the lines containing split peaks, e.g., the (hlO) line, and concentrate on 
the two peaks around the (110) position. The peak positions k a (k is the deviation 
of the wavevector from the (110) position along the (hlO) line) are determined by the 
condition da/dk = which gives 



dV fl£ 
2 = T — . 

dk dk 



(149) 



Eq. ( |149| ) means that the k-dependence of E leads to the shift 5 k = k a — ky of the 
peak position with respect to the position ky of the corresponding minimum of V(k) 
(the latter is defined by the condition dV/dk = 0). Furthermore, the right-hand side 
of Eq. (|149 ) is a function of T, while its left-hand side is T- independent. The ot(k) 
peaks will therefore change their positions with temperature. The "local" temperature 
behavior of the splitting is reflected in the sign of the derivative dxki which can be 
calculated by expanding the derivatives in Eq. ( |149| ) in powers of small changes of T 
and kj and retaining only linear terms: 



dk a 
~dT 




dkdT J 



2 &v_ T d 2 ^ 



dk 2 



dk 2 



(150) 



Our aim now is to develop a kind of minimal, i.e., simplest possible, theory which would 
be able to describe essential features of the considered effect (and, apropos, two other 
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anomalies discussed in this section). Interestingly, the approach formally rather similar 
to the Landau theory of second-order phase transitions 7 could be used. Indeed, in the 
case of not very large splitting the expansion of V(k) and in powers of k can 

be used. To describe the split minimum of V(k), only the second- and fourth-order 
terms are necessary; since the (110) position serves as the origin, the expansions do not 
contain odd powers of k. We therefore assume that in the area of the splitting V{k) 
and have the following approximate form, 

V(k) = V(0) + l -A v k 2 + ^B v k 4 , (151) 
t(k) = E(0) + ^Axk 2 + ^Bsfc 4 , (152) 

where Ay < 0, By > 0, the signs of A-% and B% are arbitrary (there are no apparent 
restrictions on the behavior of the self-energy), and k = corresponds to the (110) 
position. The resulting inverse intensity a~ l (k) has exactly the form of the Landau 
free energy functional in the low-symmetry phase where the latter possesses a double 
minimum; this implies 2|v4y| + TA-£ > and 2By — TB-% > 0. The wavevector k 
plays the role of the order parameter. Substituting approximations ( |151|) and (|152|) 
into general Eqs. ( |149| ) and ( |150| ), we get 



2\Ay\+TAy 

k « = MhTB > ( 153 ) 



, i dk a 1 

a Qj, - 2 



2B V - TBj: 
Ax + TdAy/dT Bx + TdB s /dT 
2\A V \+TA S + 2B V -TB J] 



(154) 



while ky = ±y\Ay\/By. It is seen that the shifts of the two peaks and their temper- 
ature derivatives have opposite signs and the same absolute values. Eq. (|154|) clearly 
shows two scenarios for the temperature behavior of the splitting, depending on the 
sign of its right-hand side which can be either positive or negative. The first one is 
the increase of the splitting with temperature discussed above; this takes place when 
the right-hand side of Eq. ( |154| ) is negative. Apart from that, the theory predicts that 
the decrease of the splitting with increasing temperature is also possible. This regime 
corresponds to the case of positive right-hand side of Eq. (J154|), and such temperature 



dependence has not yet been observed experimentally. Thus, the behavior of the self- 
energy determines whether the splitting increases or decreases with temperature. At 
high temperatures the correction to the wavevector- independent Y I RPA (Eq. ( |122|) ), 



and, therefore, A^ and B^, are of order A 2 (i.e., T~ 2 ). From Eq. (|153|) it then follows 
that the absolute value of the shift 8k decreases as T~ l with temperature, unless the 
corresponding coefficient identically vanishes. 



Coalescence of Fermi surface-related intensity peaks 



The analogy with the Landau theory of phase transitions, though rather formal, 
immediately leads to the following question: in the considered case, what would cor- 
respond to the transition point? The answer to this question is fairly obvious; there 
exists a possibility for the splitting to disappear at some point as temperature decreases, 
before the transition to the low-symmetry phase occurs. 44 This anomaly was neither 
observed experimentally nor correctly described theoretically, though the coalescence 
of intensity peaks (unrelated to any Fermi surface effects) was found for the exactly 
solvable ID Ising model 45 and in the CVM calculations for the 2D ANNNI model. 46 
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As we will show in this subsection, Eq. ( f75|) provides clear understanding of how such 
effect takes place. 

In the treatment given above it is, in fact, implicitly assumed that the wavevec- 
tor dependence of the interaction term 2{3V(k) in Eq. (|7^) dominates, at least in the 
area of the splitting, i.e., near the (110) position. In this case the shape of the diffuse 
intensity closely follows that of V^(k); in particular, there exists one-to-one correspon- 
dence between the split minima of the interaction and the split intensity peaks. The 
variation of the self-energy with k in this part of the reciprocal space is comparatively 
weak, though qualitatively important for the description of the temperature-dependent 
splitting. This assumption is certainly correct at sufficiently high temperatures, where 
the RPA (KCM) approximation in which the self-energy is k-independent works rea- 
sonably well. Meanwhile, as temperature starts to decrease, the correction to the RPA 
self-energy (oc T~ 2 ) grows faster than the interaction term 2(3V(\<l) (oc T _1 ). We can 
then encounter a situation when the variations of S(k) and 2f3V{\<L) with the wavevector 
are of the same order of magnitude. With temperature further decreasing, the wavevec- 
tor dependence of the self-energy can even become dominant. In this regime positions 
of the diffuse intensity peaks would be determined by the maxima of S(k). 
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Figure 2. Schematical temperature dependence 
in the coalescence case. 
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the intensity profile (a) and the peak positions (b) 



The behavior of the self-energy in the k-space is, in general, qualitatively different 
from that of the interaction. In particular, there is no special reason to expect that the 
self-energy would have any extrema away from the special points. Let us assume that 
the self-energy does not have such extrema and that the variation of the self-energy in 
the reciprocal space becomes more and more important in comparison with the corre- 
sponding variation of 2(3V{k) as temperature decreases. Then the qualitative picture 
of the temperature behavior of the splitting is as follows. At high temperatures the 
self-energy is almost k-independent, and the intensity peak positions deviate little from 
those of the corresponding minima of l^(k). As temperature decreases, the wavevec- 
tor dependence of £ becomes more pronounced; the peak positions move farther away 
from the positions of ^(k) minima and towards those special points k at which S(k) 
has maxima. Eventually, as temperature reaches certain value T , the intensity peaks 
coalesce at these special points and the splitting disappears (Figure 2). 

The coalescence temperature Tq can be found from the condition of vanishing second 



derivative a" = d 2 a/dk 2 of the intensity with respect to the wavevector at ko; the sign 
of this derivative controls the presence or absence of the splitting. At the special points 
all the first derivatives are equal to zero, and from Eq. ([75]) we obtain 



a"(k ) = c(l - c)a 2 (k ) [E"(k ) - 2f3V"(k )\ . (155) 

Therefore, the splitting disappears when the second derivatives, or curvatures, of the 
self-energy and of the interaction term 2/3V(k) at the special point ko compensate each 
other, i.e., when 

E"(k ) = 2f3V"(k ) . (156) 

To analyse the behavior of the splitting close to the coalescence point, it is convenient 
to use the same Landau-type approach as in the previous case. In this temperature range 
the splitting above T is small, and expansions (|151|) and ( |152| ) are valid. Substituting 
them into Eq. (|75|), we get 

a -\k) = a~\0) + ^Ak 2 + ^Bk i , (157) 
where second-order coefficient 



A = c{l-c)(-Av + 2/3A v ) (158) 



vanishes at T = T (see Eq. ( p.56| )), while the fourth-order coefficient B remains positive 



at that temperature. We can then, at temperatures close to T , regard A as linear in 
T — T with a negative coefficient and B as temperature- independent. Thus, the inverse 
intensity a~ 1 {k) behaves almost in the same way that the Landau free energy. The only 
difference here is that the role of temperature is reversed; a~ l (k) has a double minimum 
above the coalescence temperature To and a single minimum below it. Therefore, at 
small positive values of T — T the splitting increases with temperature as (T — To) 1 / 2 . 
Contrary to the corresponding result of the genuine Landau theory, obtained bifurcation 
exponent is exact, since the intensity is an analytical function of the wavevector and 
can legitimately be expanded into the Taylor series. At higher temperatures behavior 
of the splitting changes, and sufficiently far away from T it starts to approach the value 
dictated by the interaction V(k). 



"Thermal" splitting of intensity peaks 



We have just considered the situation when the double-well profile of the interaction 
in the vicinity of a special point is compensated by the wavevector dependence of the 
self-energy which has a simple maximum at this position. The competing curvatures 
2j3V"(k ) and S"(fc ) are both negative. As a result, the second derivative (|155|) of the 
intensity vanishes at some temperature, and the splitting induced by the Fermi surface 
effects disappears. 

It is not very difficult to realize that another kind of curvature compensation is pos- 
sible; this is the case in which both curvatures are positive. 39 This is, perhaps, the most 
physically interesting situation: here the interaction with a single minimum produces 
the intensity peak with no fine structure at higher temperatures (in full accordance with 
the RPA-like considerations), but then, as temperature decreases and the wavevector 
dependence of the self-energy becomes more and more significant, the compensation 
takes place and the intensity peak splits. This is especially probable when the min- 
imum of V(k) is shallow (i.e., V"(k ) is small); in particular, in the limiting case of 
vanishing V"(ko) it is the curvature of that controls the fine structure (single- vs. 



no 



double-peak) of the maximum of a(k). The application of the Landau-type description 
gives essentially the same results as before. However, the coefficient in front of T — T in 
A is now positive; the inverse intensity a~ 1 (k) has a single minimum above the splitting 
temperature T and a double minimum below this point, and the splitting increases as 
(T — T) 1 / 2 with decreasing temperature at small negative T — T values (Figure 3) . 
The bifurcation exponent is again exact, for the reasons mentioned above. 




Figure 3. Schematical temperature dependence of the intensity profile (a) and the peak positions (b) 
in the case of the "thermal" splitting. 



This type of behavior was recently experimentally observed (and subsequently re- 
produced in the MC simulations) for the Pt-V alloy system by Le Bolloc'h et al. 37 In 
this system the splitting of the (100) intensity peak along the (hOO) line occured with 
decreasing concentration of vanadium rather then temperature. The explanation of 
this anomaly was proposed by the author, 39 who also pointed out that similar splitting 
should take place when temperature decreases at fixed composition. The predicted 
effect was then discovered in the MC simulations by the same group. 47 Experimental 
confirmation of its existence remains to be seen. 
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